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[l] We present the first observations of large amplitude waves 
in a well-defined electron diffusion region based on the criteria 
described by Scudder et al.[ 2012] at the subsolar magnetopause 
using data from one Time History of Events and Macroscale 
Interactions during Substorms (THEMIS) satellite. These 
waves identified as whistler mode waves, electrostatic solitary 
waves, lower hybrid waves, and electrostatic electron 
cyclotron waves, are observed in the same 12 s waveform 
capture and in association with signatures of active magnetic 
reconnection. The large amplitude waves in the electron 
diffusion region are coincident with abrupt increases in 
electron parallel temperature suggesting strong wave heating. 
The whistler mode waves, which are at the electron scale 
and which enable us to probe electron dynamics in the 
diffusion region were analyzed in detail. The energetic 
electrons (~30keV) within the electron diffusion region 
have anisotropic distributions with T ei /T e || > 1 that 
may provide the free energy for the whistler mode waves. 
The energetic anisotropic electrons may be produced 
during the reconnection process. The whistler mode waves 
propagate away from the center of the “X-line” along 
magnetic field lines, suggesting that the electron diffusion 
region is a possible source region of the whistler mode waves. 
Citation: Tang, X., C. Cattell, J. Dombeck, L. Dai, L. B. Wilson III, 
A. Breneman, and A. Hupach (2013), THEMIS observations of the 
magnetopause electron diffusion region: Large amplitude waves and 
heated electrons, Geophvs. Res. Lett., 40, 2884—2890, doi:10.1002/ 
grl.50565. 

1. Introduction 

[ 2 ] Magnetic reconnection is considered to be an important 
energy conversion process that occurs in a variety of plasma 
environments. At the Earth’s magnetopause, it facilitates the 
entry of solar wind plasma and electromagnetic energy into 
the magnetosphere. Strong wave turbulence covering a broad 
range of frequencies has been frequently observed near the 
reconnection site [Retind et al., 2006; Khotyaintsev et al., 
2006; Cattell et al., 2005]. More specifically, these wave 
modes include the whistler mode (WH) waves [Deng and 
Matsumoto, 2001; Petkaki et al., 2006], electrostatic solitary 
waves (ESWs) [Farrell etal., 2002; Matsumoto etal., 2003], 
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lower hybrid (LH) waves [Cattell and Mozer, 1986; Bale 
et al., 2002], kinetic Alfven waves [Chaston et al., 2005], 
and Langmuir/upper hybrid waves [Farrell et al., 2002], 
The effect of different wave modes on the reconnection 
process has been a problem of long-standing interest — for 
their role in anomalous resistivity, particle acceleration, 
energy transport, and formation of reconnection sites [Huba 
et al., 1977; Labelle and Treumann , 1988; Treumann et al., 
1991; Drake et al., 2003; Roytershteyn et al., 2012], 

[ 3 ] The physics of reconnection occurs on the larger spatial 
scale (ion diffusion region) and on the smaller scale (electron 
diffusion region) that are associated with the ion and electron 
dynamics, respectively. Properties of collisionless asymmet- 
ric reconnection have been compared to those of symmetric 
reconnection on ion-scale physics [Pritchett, 2008] and elec- 
tron-scale physics [Mozer and Pritchett, 2011, and references 
therein]. Simulations have been performed to investigate the 
effects of diamagnetic drift [Swisdak et al., 2003], guide 
magnetic field [Pritchett and Mozer, 2009], and dissipation 
mechanism [Aunai et al., 2013] on collisionless asymmetric 
reconnection. The efficiency of collisionless asymmetric 
reconnection is controlled by the physics in the electron dif- 
fusion region (EDR). Observations of EDRs have been made 
by Polar at the subsolar magnetopause [Mozer et al., 2003], 
by Wind in the magnetotail [Qieroset et al., 2001], and by 
Cluster in the magnetosheath [Phan et al., 2007]. Recent 
simulations and observations of EDRs during collisionless 
antiparallel reconnection in Earth’s magnetotail [Xg et al., 
2011] report that the diffusion region is characterized by a 
narrow extended layer containing electron jets. It is shown 
that the jets in the layer are driven by electron pressure anisot- 
ropy P e || » P eJ _, and the anisotropy is responsible for the 
structure of the EDR [Xg et al., 2011], Mozer [2005] has 
identified EDRs on the basis of the nonzero parallel electric 
field, a large perpendicular electric field compared to the 
reconnection electric field, a large electromagnetic energy 
conversion rate and accelerated electrons, and a topological 
boundary that separates regions having different E x B IB 2 
flows with thickness of the order of the electron skin depth. 
In this paper, we concentrate on a specific way to identify 
the EDR described by Scudder et al. [2012] who report 
spatially resolved diagnostic signatures of a demagnetized 
EDR observed by Polar at the Earth’s magnetopause. The 
criteria for being in an EDR should be electron specific since 
electrons are locally disrupted in the EDR. The four dimen- 
sionless scalar diagnostics that were used to find the EDR 
are peak electron thermal Mach numbers M ei = |U f |/(w fi ) 
>1.5 (where U e represents electron bulk velocity and {w eA _) 
is the electron thermal speed derived from the average 
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perpendicular temperatures), electron temperature anisotropy 
An e = T e || l(T e ± ) > 7, calibrated agyrotropy of electron pres- 
sure tensor A0 e =2|l — a|/(l+a)> 1 (wherea^Pe^j/Pej,^), 
and strong ( ^ 150 eV) increases in electron temperature 
[Scudder et al., 2012]. The electron pressure agyrotropy is 
a measure of the departure of the pressure tensor from 
cylindrical symmetry about the local magnetic field and it 
reflects the demagnetization of the thermal electrons, which 
enables collisionless reconnection. 

[ 4 ] This paper focuses on the WH waves, which are an 
important candidate for the anomalous resistivity, particle 
acceleration, and heating. WH waves may be driven unstable 
by superthermal electrons with temperature anisotropies of 
T e ±/T e || > 1 in the magnetosphere [Kennel and Petschek , 
1966] and current-driven plasma instabilities [Gurnett 
et al., 1976] or energetic electron beams [Zhang et al., 
1999] in the magnetotail. WH waves are one of the most 
ubiquitous wave modes observed in space plasmas. 
Observations of WH waves at the Earth’s magnetopause 
have been made by Deng and Matsumoto [2001]. Electron 
anisotropy, due to compression of the magnetopause or LH 
drift waves, may be the generation mechanism of WH waves 
in the magnetopause current sheet [Karimabadi et al., 2004], 
WH waves in the EDR may play a significant role in the 
microphysics of reconnection as they are excited on electron 
scales. It is believed that WH waves in the magnetopause 
current sheet may affect the instability of the current 
sheet to reconnection via tearing. The generation of the out- 
of-plane component of the magnetic field is suggested to 
be a signature of whistler-mediated reconnection [Mandt 
et al., 1994]. It has also been suggested that the strongest 
whistler emissions are observed on the most recently opened 
magneto spheric flux tubes due to magnetic reconnection 
[ Vaivads etal., 2007]. One recent simulation study concludes 
that WH waves do not control the dissipation processes 
of reconnection but are generated as a result of the 
reconnection processes [Fujimoto and Sydora, 2008], In this 
paper, we present an example of a reconnection event at the 
subsolar magnetopause observed by Time History of 
Events and Macroscale Interactions during Substorms 
(THEMIS). In section 2, we describe the data sets and 
analysis techniques. In section 3, we show the observations. 
Finally, we discuss the conclusions of our study in 
section 4. 

2. Data Sets and Analysis 

[ 5 ] The THEMIS mission consists of five identically 
instrumented spacecraft. The Electric Field Instrument (EFI) 
measures three components of the electric field [Bonnell 
et al., 2008]. The instrument provides continuous coverage at 
128 samples/s in survey mode and waveform captures at 
8192 samples/s in wave burst mode. The Magnetic Fields 
Experiment (MFE) includes a fluxgate magnetometer (FGM), 
which measures DC magnetic field with a sampling rate of 
128 samples/s in the high rate mode or 4 samples/s in the low 
rate mode, and a search coil magnetometer, (SCM) which mea- 
sures magnetic fluctuations sampled at 8192 samples/s in the 
burst mode [Roux et al., 2008]. Particle data are measured by 
the Electrostatic Analyzer (ESA) [McFadden et al., 2008] 
and by the Solid State Telescope (SST) [Angelopoulos, 
2008]. The ESA measures plasma over the energy range of a 
few eV up to 30 keV for electrons and 25 keV for ions. The 


SST measures the distribution functions of superthermal 
particles within the energy range from 25 keV to 6 MeV for 
electrons and 900 keV for ions. Because the measured quasi- 
static electric field component along the spin axis has large un- 
certainty due to the short boom along the spin axis, we use 
EB = 0 to determine the electric field used to calculate the 
E x B IB 2 velocity. The coordinate systems used in this paper 
include geocentric solar magnetospheric (GSM) coordi- 
nates and field-aligned coordinates (FAC). The FAC is 
defined in the following way: The positive Z axis points in 
the direction of the magnetic field at the spacecraft’s 
location. The positive X axis lies in the plane of the 
magnetic field line passing through the spacecraft’s location, 
perpendicular to the Z axis, and points inward (toward the 
inside of the field line). The positive Y axis completes the 
orthogonal right-handed system. Waveforms are analyzed 
dynamically in time and frequency using Morlet wavelet 
transform. The wave vector is determined using Minimum 
Variance Analysis (MVA) on band pass-filtered three 
magnetic field components. 

3. Observations 

[6] Figure 1 shows a 5.5 min interval of the field and 
plasma observations made by probe E of the THEMIS 
mission on 27 August 2009. The boundary normal direction 
(determined from MVA on the ambient magnetic field) was 
[0.99, 0.015, —0.12] in GSM coordinates and almost 
identical to the GSM-X direction, consistent with the 
spacecraft being near the subsolar point (indicated by the 
position parameters at the bottom of Figure 1). The 
spacecraft travels from the outer magnetosphere (SP) 
through the magnetopause (MP, indicated by two light 
green shaded bands) into the magnetosheath (SH). The 
purple shaded band shows a ~12 s interval of magnetic 
burst data capture. 

[ 7 ] The magnetopause crossing can be seen in the change 
in Bz from positive to negative (Figure la). The differential 
energy flux of electrons in Figure lj shows that in the magne- 
tosphere, where Bz was positive, higher-energy electrons 
were encountered; while in the magnetosheath, where Bz 
was negative, lower-energy electrons were measured. The 
spacecraft passed from the lower plasma density magneto- 
sphere to the high-density (factor of 100) magnetosheath 
via a region of mixed magnetosheath/magnetospheric 
plasmas comprising the low-latitude boundary layer as 
shown in Figure lg. This observation is evidence for the 
transport of solar wind plasma across the magnetopause. 
The presence of accelerated plasma flow is seen through the 
magnetopause current sheet as shown in Figures Id and le. 
The spacecraft crossed the magnetopause south of the separa- 
tor, as suggested by the negative GSM-Z component of ion 
flow velocity and the result of Walen test [Sonnerup et al., 
1981], Based on the prediction of the magnetopause 
reconnection model being a rotational discontinuity 
[Paschmann et al., 1979], the Walen test states that the 
observed flow velocity change between a point in the magne- 
topause and a reference point in the adjacent magnetosheath 
equals the predicted modified Alfven velocity change. The 
ratios of magnitudes between the observed flow velocity 
change and the theoretically predicted modified Alfven 
velocity change range from 0.30 to 0.66, which is reasonable 
within the range of uncertainty. The angle deviations 
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Figure 1. A reconnection event at the subsolar magnetopause observed by THEMIS-E on 27 August 2009. (a) 4 samples/s 
magnetic field data in GSM. (b) Perpendicular X component of the burst magnetic field at 8192 samples/s in FAC. 
(c) Perpendicular Y component of the electric field at 128 samples/s in FAC. (d and e) Ion and electron bulk flow velocity 
in GSM, respectively, (f) Comparisons of GSM-Z component of the E x BIB 2 velocity with the GSM-Z component of ion 
and electron perpendicular flow velocity with respect to the ambient magnetic field, respectively, (g) Ion and electron densi- 
ties. (h) Electron temperatures, (i) Electron agyrotropy and Mach number (scale to the left) and temperature anisotropy (scale 
to the right), (j) Differential energy flux for electrons measured by ESA. (k and 1) Electron pitch angle spectra for lower-energy 
(7 eV-26keV) electrons measured by ESA and higher-energy (3 1 719 keV) electrons measured by SST, respectively. 


between the observed and the predicted velocity changes are 
almost 180° for this event. This antiparallel relation indi- 
cates that the normal magnetic field component Bn should 
be positive and thus the spacecraft crossed south of the 
separator [ Sonnerup et ai, 1981]. This result is consistent 
with the positive polarity of the GSM-X component of the 
magnetic field, which is an approximation for Bn shown 
in Figure la (around 15:35:30 UT). The encounter with 
the magnetopause current sheet is associated with fast ion 
jetting consistent with the Walen relation and fast electron 
flows, indicating that reconnection is occurring. Magnetic 
reconnection is generally considered to be the primary 
mechanism through which transport of plasma and energy 
across the magnetopause occurs. 


3.1. Identification of Electron Diffusion Region 

[8] Enhanced wave activity can be seen associated with the 
magnetopause crossing from Figures la (ambient magnetic 
field), lb (burst magnetic field), and lc (electric fluctuations). 
We note that the electric fluctuations maximize during the 
magnetic burst interval. During this interval, electron flow 
speed (Figure le), anisotropy, agyrotropy, and Mach number 
(Figure li) also maximize. These enhanced amplitudes are coin- 
cident with abmpt increases in electron parallel temperature T e y 
shown in Figure lh suggesting strong wave heating. As will be 
discussed in more detail in section 3.3, the observed intense 
waves may provide the observed electron heating. All these fea- 
tures, along with the fact that the electron perpendicular flow 
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Figure 2. Example of WH waves within the EDR. (a, b, and c) A 0.53 s interval in FAC of perpendicular A component of the 
burst magnetic field, associated wavelet power spectrum, and whistler Poynting flux, respectively, (d and e) Expanded views 
of the filtered whistler waveforms of the perpendicular X component of the burst magnetic field and the perpendicular Y com- 
ponent of the electric field over the time interval indicated by the purple bar in Figure 2a. (f) Hodogram of the filtered burst 
magnetic field waveforms in FAC for the interval indicated by the green bar in Figure 2d. The black star and dot mark the 
beginning and ending of the wave field, respectively, (g and h) Distribution functions of lower-energy electrons (up to 
20keV) measured by ESA and higher-energy electrons (20-700 keV) measured by SST observed at times close to the WH 
waves, respectively. The horizontal axis is parallel to the ambient magnetic field and the bulk velocity defines the plane. 


velocity is not consistent with the E x B IB 2 velocity during the 
magnetic burst interval (Figure If), provide evidence for the 
detection of an EDR. These signatures are consistent with 
the simulation and observations of Scudder et al. [2012]. It 
should be noted that the small peaks in the electron thermal 
parameters correspond to entries of topological boundaries of 
the magnetic field. For example, the peaks around 15:33:30 


UT are due to the entry of magnetopause boundary layer. A 
density cavity can be seen in Figure lg (purple shaded band). 
The density inside the density cavity decreases to 16% of the 
magnetopause boundary level (30 cm -3 ), which might indicate 
the center of the EDR. This density cavity is not a signature of 
magnetospheric separatrix as the density inside the density 
cavity does not drop with respect to its magnetospheric level. 
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Figure 3. Example of electrostatic waves within the EDR. (a, b, and c) A 0. 1 s interval of the perpendicular A component of 
the burst magnetic field, the parallel Z component of the electric field waveform capture and the associated wavelet power 
spectrum of the parallel electric field in FAC, respectively, (d, e, and f) Flodograms of the electric field waveforms with an 
interval indicated by the light blue arrows below Figure 3c. (g) Fourier wave power versus frequency with the same time 
interval as the hodograms. 


From the deviations of ion perpendicular flow velocity from 
the E x B/B 2 velocity (Figure If), it can be seen that the 
ion diffusion region was encountered probably during the 
following two time intervals: 15:34:40-15:36:20 UT and 
15:36:50-15:37:10 UT. 

3.2. Distributions of Electrons 

[ 9 ] Figures lk and 11 show the electron pitch angle spectra 
for lower-energy electrons measured by ESA and higher- 
energy electrons measured by SST, respectively. The 
magnetic field (Figure la) and electric field (Figure lc) fluc- 
tuations enhance in the magnetopause boundary layer and in 
the magnetosheath with field-aligned and counter-streaming 
lower-energy electrons as shown in Figure lk. However, 
the higher-energy electron pitch angle (Figure 11) enhances 
around 90° in the boundary layer and in the magneto sheath, 
especially during the purple shaded magnetic burst interval. 
Distinct from the electron distributions in the boundary layer 
and in the magnetosheath, the higher-energy electron pitch 
angle enhances at 0° and 180°, and the lower-energy elec- 
trons are more isotropic near the current sheet center around 
15:36:00 UT and 15:37:00 UT. The lower-energy electrons 


in the identified EDR in Figure lk can be seen streaming 
away from the “X-line” (antiparallel to the local magnetic 
field). These electrons are associated with Hall current paral- 
lel to the magnetic field. This indicates that the identified 
EDR is close to the “X-line.” 

3.3. Observations of Waves 

[ 10 ] Figure 2 shows an example of the identified WH waves 
at the time indicated by a black vertical line in Figure lb. It can 
be seen from Figures 2a and 2b that the waves have frequen- 
cies from 0.1 to 0.6 f ce (electron cyclotron frequency) with 
amplitudes up to 3 nT (peak-peak). As can be seen from 
Figure 2c, the wave Poynting flux is mostly antiparallel to 
the ambient magnetic field. The wave vector (£=[— 0.23, 
—0.09, 0.97] in FAC), determined from MVA on band pass- 
filtered wave magnetic field, is nearly along the background 
magnetic field. The wave propagation angle with respect to 
the ambient magnetic field 0 kH is determined as -166° since 
the Poynting flux is mostly antiparallel to the magnetic field. 
An expanded view of the WH waves can be seen in 
Figures 2d and 2e, which respectively show the filtered 
(200-2000 Hz) waveforms of the burst magnetic and electric 
field data over the time interval indicated by the purple bar 
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in Figure 2a. Figure 2f shows that the WFI waves are circularly 
right-handed polarized with respect to the ambient magnetic 
field as expected. The electron distribution functions observed 
at times close to and/or concurrently with the WFI waves are 
shown in Figures 2g and 2h. The lower-energy electrons 
(~100eV) shown in Figure 2g have anisotropic distributions 
with T e JT e \\ < 1. Flowever, the energetic electrons 
(~30keV) shown in Figure 2h have anisotropic distributions 
with a larger population moving in perpendicular direction. 
Broad-banded emissions with strong electric (~10mV/m) 
and magnetic (—40 nT) field fluctuations below ion cyclotron 
frequency (~1 FIz), and electric (~30mV/m) and magnetic 
(~20 nT) field fluctuations below LIT frequency (~30Hz) are 
also detected during this magnetopause crossing (not shown). 
These intense wave emissions may provide the observed elec- 
tron heating associated with the magnetopause crossing. 

[ 1 1 ] F igure 3 shows examples of electrostatic waves at a time 
preceding the waves in Figure 2 by 0.05 s. WH waves at a 
frequency of ~0.3 f/ e are observed in the magnetic fluctuations 
in Figure 3a from 15:35:33.650 UT to 15:35:33.700 UT. This 
time interval overlaps with that of the ESWs, suggesting a 
possible coupling of WFI waves and ESWs. ESWs (up to 
30mV/m) indicated by the magenta arrows in Figure 3b have 
a broad spectrum, which extends from 200 to 3000 FIz shown 
in Figure 3c. The high-frequency electrostatic waves (up to 
35 mV/m) labeled by light blue arrows in Figure 3b have 
power that peaks at f ce , which can be seen from both the wave- 
let power spectrum in Figure 3c and the Fourier power in 
Figure 3g. During this time interval, there is no power in the 
wave magnetic field at f ce (not shown). This wave mode is 
linearly polarized, as shown in the hodograms in Figures 3d, 
3e, and 3f in FAC with an interval indicated by the light blue 
arrows below Figure 3c. Occasionally, these high-frequency 
emissions are seen associated with clear harmonics, possibly 
suggesting electrostatic electron cyclotron waves. 

4. Discussion and Conclusions 

[ 12 ] We have presented the first observations of intense 
waves in the EDR in a subsolar magnetopause reconnection 
region. The identification of the EDR in this event is based 
on the occurrence of signatures of strong electron heating, 
large electron thermal anisotropy, agyrotropy, and Mach 
number and electron velocity not consistent with the E x B 
velocity, consistent with the simulation and observations of 
Scudder et al. [2012]. The lower-energy electrons 
(-100 eV) with anisotropic distributions of T e JT e \ i \ i < 1 
within the EDR may have been heated by the observed waves 
with frequency below the LFI frequency, consistent with the 
suggestion that LH waves lead to electron heating in the par- 
allel direction [Cairns and McMillan, 2005], 

[ 13 ] We identified intense WFI waves inside the EDR. This 
is inconsistent with reported simulation results that indicated 
WH waves are only driven downstream of an EDR [Fujimoto 
and Sydora, 2008]. The WH waves seen by THEMIS in the 
EDR propagate almost antiparallel to the ambient magnetic 
field, and the Poynting flux indicates that the WH waves 
propagate away from the center of the “X-line” along mag- 
netic field lines. The observed electron temperature anisot- 
ropy of T e j_/T e || > 1 for energies above 20 keV may be the 
source of free energy for the generation of the WH waves. 
The energetic electron anisotropy may be produced by adia- 
batic heating in the perpendicular direction as the locally 


intensified magnetic field can accelerate electrons in the per- 
pendicular direction [Fujimoto and Sydora, 2008]. On the 
field lines directly connected to the EDR, the energetic elec- 
tron anisotropy may also be due to higher-energy field- 
aligned electrons (accelerated by the reconnection process) 
being lost to the magnetosheath [Stenberg et al., 2005], 
WH waves can scatter the electrons in pitch-angle distribu- 
tion and relax the temperature anisotropy. Studies of large 
amplitude whistlers in the inner magnetosphere have pro- 
vided evidence for rapid scattering and/or energization 
[Cattell et al., 2008]. WH waves may play a significant role 
in the microphysics of reconnection through the enabling of 
a current sheet instability, the decoupling of electrons, the ac- 
celeration and heating of particles, and the transport of en- 
ergy away from the reconnection region. 

[ 14 ] A possible coupling of electrostatic electron cyclotron 
waves and ESWs with WH waves is often seen during magne- 
topause reconnection. The growth of the electrostatic waves 
may reduce the electron temperature anisotropy and reduce 
the growth rate of WH waves. The physics of wave coupling 
process is important to understand the effect of wave-wave 
interactions on the reconnection process and will be investi- 
gated in a future study. This study provides further evidence that 
the plasma waves can play a significant role in the microphysics 
of magnetic reconnection at the Earth’s magnetopause. 
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